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Abstract 

We have found some new exact static spherically symmetric interior solutions of metric f{R) gravi- 

o 

. tational theories describing the equilibrium configuration of a star. Then the solution is matched to the 

^ 'l exterior solution and thus gives a complete description of a star in R^~"^^ theory. 
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1 Introduction 

One of the important problems in the gravitational theories is finding solutions under the assumption of 
spatial isotropy and time independence, i.e. static spherically symmetric solutions. Then the results can 



^ ■ be used to treat the solar system tests and also for computing the pressure, density and gravitational field 

o : 

0^ ■ within a spherically symmetric static star. Regarding the standard general relativity, there is a large number 



of exact solutions of Einstein's field equations for perfect spherically symmetric static fluids in the literatures. 
For a complete list see [I]. From these solutions which are in closed algebraic form, only eleven solutions 
satisfy physical criteria mentioned in |T and therefore are physically acceptable, (nine solutions in fl], one 
^ . in [2] and one in [3]). 

A particular class of gravitational theories which has been received great interest recently, is /(i?) modified 
theories of gravity [4]. In [5j spherically symmetric vacuum solutions of f{R) theories are studied. In particular 
it is shown that the modified gravitational theories have the general static blackhole solutions with non- 
vanishing cosmological constant, Schwarzschild-deSitter solutionlSl [6l [?]• Considering scalar-tensor theories 
of gravity in the Einstein frame, a reconstruction method is introduced in ^ to generate some spherically 
symmetric vacuum space-times. Since f{R) gravity is dynamically equivalent to Brans-Dicke theory with zero 
coupling constant, this method is applicable to f{R) gravity in the Einstein frame. Also in the recent years 
some exact solutions for static spherically symmetric perfect fluid in f{R) theory of gravity are presented 
[S1[TUJ[TT]. Since the field equations of f{R) theory are differential equations of fourth order, more integration 



constants need to be determined in order to get a unique solution. This shows that in general, the exterior 
space-time geometry is not described only by the mass of the star, and thus the validity of BirkhofF theorem 
is lost[TH[T3]- Other details of mass distribution, i.e. the gravitational multipole moments, are needed to 
obtain a unique solution for both interior and exterior regions of a star. As a result the exterior solutions of 
f{R) theory are not necessarily asymptotically flat. 

Another issue related to the star evolution is the singularity problem in f{R) gravity which has been first 
introduced in |14] and then disscused for stellar configurations in |15l 116] which maybe treated by adding 
higher powers of R [17] . The classification of future singularities according to their properties into four classes 
is studied in [TH]. These singularities occur in scalar or fiuid models of dark energy and in /(i?) or /(G) 
(G denotes the Gauss- Bonnet invariant) [19] or other models leading to an accelerating universe. The fact 
that all such singularities is not unique for f{R) but occur in all known dark energy models is explained in 

Moreover the consistency condition of any f{R) theory with the solar system tests is the occurrence 
of the chameleon mechanism [^H [H]. It is one way to avoid the Doglov-Kawaski instability [53^. Also a 
positive value of F{R) — df{R) /dR is needed to guarantee the positivity of the effective gravitational constant 
Ge// = G/F{R) U. 

In [121 [13], by considering the Lagrangian of the form of R^^^ , two exact spherically symmetric vacuum 
solutions are presented which do not have the asymptotic behavior of the Minkowski space-time. It seems 
that the non-existence of asymptotic flatness is one of the reasons to study weather /(i?) theory can explain 
the flat rotation curve of galaxies or not [HI [Ml HSHH] • Moreover f{R) static spherically symmetric solution 
for stellar systems has been investigated in [27]. The authors show that specifying the pressure and the 
density of matter do not uniquely determine the functional form of f{R)- As mentioned before, this is 
because of appearing higher order derivatives of quantities in the equations which require matching of higher 
order derivatives of variables at the surface of star. It is shown that assuming Schwartzschild-deSitter metric 
outside the mass distribution as a boundary condition, induces tight constraints on the quantities [57]. 

In this paper, we discuss the gravitational equations of f{R) theory for the interior of a spherically 
symmetric matter distribution and derive a few new solutions in some closed algebraic form. For these 
solutions, the dependence of f{R) function to the scalar curvature can be found. We divide these solutions 
into two groups. After analyzing each metric, we see that one of the obtained metrics is appropriate for 
describing a star. At the end, we consider this solution in more detail and match it to the corresponding 
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exterior metric solution. 



2 Spherically symmetric static solutions of f{R) gravity 

The dynaniics of f{R) theory can be determined by the following action: 

S = J d^xV^{f{R)) + i^C„,) (1) 

In the metric approach, the gravitational field equations can be obtained via the variation of action with 
respect to metric: 

F{R)R^,. - lfiR)g^.. - V^V,F(i?) + (?^.V"V„F(i?) = kT^, (2) 

where T^i^ is the energy-momentum tensor of matter and F{R) = df /dR. Considering a static and spherically 
symmetric distribution of matter, the line element has the following form: 

ds^ = A{r)dt'^ - B{r)dr'^ - r'^dU'^ (3) 

and assuming the matter to be described by a perfect fluid, we have: 

T'''' = (p + p)u^'u'' - pg''" (4) 

where p{r) and p{r) are the proper mass density and the isotropic pressure of the fluid in the instantaneous 
rest frame of the fluid and is its four velocity. Substituting ([3]) and in ([5]), one obtains the static 
spherically symmetric field equations of f{R) theory with matter. Assuming a specific form of /(-R), these 
equations would be three fourth order nonlinear differential equations. Supplementing by an equation of state 
for matter: 

p = p{p) (5) 

one has a complete set of equations for four unknown quantities. A, B, p and p. To carry out the calculation, it 
is desirable to re-express the equations in a somewhat different form [S] . It can be shown that the expression: 

K[iA = (6) 

with no summation on repeated indices, is independent of the index p,. From K\^i^ — iiTj^] = 0, we have: 

2F^+rF'^-2rF" + 2rK^{p + p) = Q (7) 

where X{r) — A{r)B{r) and a prime denotes the derivative with respect to r. Also K^t-^ — Kyg^ = results in: 

f F' X'\ / 2A\ 2A 2X X , , 



Furthermore, froni the conservation equation of energy momentum tensor, one obtains: 



1 



P + P 



2p' 



(9) 



And also the corresponding curvature scalar is given by: 



i? — TT 




) 



+ 



2A 



(10) 



Equations (O, d?]), ([U, dU and ([T0| completely determine the unknown quantities A, X, p, p and R. These 
are the relativistic equations of stellar structure in f{R) theory of gravity. 

From a mathematical point of view, knowing the functional form of f{R), equation (1101) shows that the 
functions F{r) and F'{r) depend on A{r) and B{r) in a complicated and non-linear manner. Therefore we 
can not expect to obtain the analytical solutions of equations ([7]), (|S]) and © in a closed form. Although 
one could use the numerical [T^ or approximate methods such as the Newtonian approximation [10 or 
perturbation method [11 . In the Newtonian limit, assuming a polytropic equation of state of matter, p ^ , 
these equations lead to a modified version of Poisson and Lane-Emden equations |28|. 

Here we are interested in obtaining some analytical solutions in a different way. In order to do this, 
we follow Tolman's method [53] which is developed to provide a variety of explicit solutions of Einstein's 
equations for a static sphere of fluid. Assuming matter is represented by a perfect fluid, instead of using the 
equation of the state of matter, Tolman introduces a relation for the metric components or their dependence 
to radial coordinate, in order to solve the set of equations easily. Of course, it must be noted that this 
procedure leads to some expression for the density and the pressure of the perfect fluid which may not be 
physically acceptable for a star |29) . 

For /(i?) equations, it is necessary to consider two assumptions. One is the equation of the state of the 
fluid and the other is the form of f{R). Logically these two assumptions should be given to have a completely 
defined model. We should do this in another way. Two relations for X and F are proposed and then the 
equations are solved. As a byproduct of the solution, the equation of state and /(i?) are obtained. However 
it has to be noted here that not all solutions lead to physically acceptable equation of state and also that not 
for all solutions it is practical to obtain the form of f{R) analytically. 

In this paper we discuss two kinds of assumptions: X — ATq, F — For"^ (class I) and X = ATor™, F — For" 
(class II) where Xq, Fq, n and m are arbitrary constants. For any choice, a number of interesting solutions 
for A, p and p along with the corresponding forms of f{R) are obtained. To specify the integration constants, 
one must specify the boundary conditions. For a star in equilibrium, it is physical to assume that the pressure 
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drops from its central value to zero at the surface of star, Rg. Moreover the metric needs to be matched with 
the outside vacuum solution at the surface of the star, according to the Israel junction conditions. 



2.1 Class I solutions: X = X^, F = For'' 

A look at equations ([T])-®, suggests that choosing X = Xq = Constant simplifies the equations. Then on 
choosing F = For" and combining equations ([7]) and © we have 

nA' - 4n - 2 , 2Xq , , 

+ + 5 A+^=0 (11) 



with the general solution 



where 



and 



A-Ao+Air™ (12) 



— zn — 1 



= i (^1 - n ± + 14?^ - 77i2^ (14) 



and A\ is a constant. 

Finally equations ([7]) and © leads to 



p = Po+Pir™+"-' (15) 



where 



and po is a constant, and also 



where 



and 



2KXo(n + m- 2) 

In addition one can simply see that 



_ FQmn{n - l)Ai 

~ 2KXo{n + m - 2) ^ ^ 



P^-Po + Pir + P2r ^ (17) 



_ Fonjn - 1) 

- "^(n2-2n-l) ^^^^ 



Fon(n- l)(2n + 3m-4)Ai 
P2 = o..^ , ^ — (19) 



R=^ + i?2r'"-' (20) 



5 



with 



and 



— 2?i — 1 



. ("-' + ^:" + ^)-^^ (22) 

Although these are exact solutions of the /(i?) equations, but the form of /(i?) only can be reconstructed 
for special choices of n. These are listed in Table ([1]). 

Solutions I, II, III, IV, and V are not physically acceptable, because they lead to a negative sound velocity 
squared = dp/ dp. The last solution is acceptable for a thick shell. The minimum and maximum radii of 
the shell would be derived from the conditions that < < 1 and p > and p > 0. As a result we have a 
thick shell: 



3Ai . .pI 
-IT- < ^ < ^ 



<r<^ (23) 



It has to be noted here that to have a physically complete solution, this thick shell solution have to be 
matched to the vacuum solution at r = 3Ai/Ao and r ~ Pi/Pq- The Israel junction conditions for /(i?) 
theory can be found in the literature [301 Ell E2- Choosing the coordinates so that the metric is in the form: 

ds^ = N'^dz^ + hijdx'dx' (24) 

where z is in the jump direction, the Israel junction conditions are: 

A/ijj- = (25) 

Ai? = (26) 

Ai^T = (27) 

F{R)AK,j = S,, (28) 
dF{R) 



dR 



in which Kij is the extrinsic curvature defined as: 



-Aa^i? = NS (29) 



K is its trace, Kij is the traceless part of Kij, Sij is the energy-momentum tensor on the boundary (the thin 
shell), S is its trace and Sij is the traceless part of Sij. 
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Table 1: The first class of solutions with X = Xq and F — For" 



A{r) 



p{r) 



p{r) 



R{r) 



f{R) 



^s = tp 



acceptable? 



I (n = 0, m = -1) 

II (n = 0, m = 2) 
III (n = 1, TO = -2) 

IV (n = 1, TO = 2) 

V (n = 2, TO = 1) 
VI (n = 2.5, TO = -1) 



2 ^ 



4^ + Air 



constant 
constant 
constant 
constant 



_ FpAir 
PO kXo 



15FoAi 



constant 
constant 
constant 
constant 



3FoAir 

K,Xo 



^° re 2reXo->/F 







1 



12Ai 



1 

7^ 



6Ai 
Xor 

10 



any regular 
any 
fo - 2FoV-R 
fo — 2FoY^ 

fo 

/o- 



l l2Ai 
Xo 



-R 



36^1 Fq 

XnR 



10Fo(-10fl)~^^'' 
4 



-Ai 



no 



no 



no 



no 



no 



2Ai+4Xor 



thick shell (see Discussion section) 



2.2 Class II solutions: X = Xgr™, F = For"" 
Choosing 

F = For" (31) 

and 

X = Xor™ (32) 

the equations dT!)-© can be solved leading to: 

A = Aor"' + + A2r^' (33) 

where 

Ao - Xo5; 5 - — — ^ ~\ ^ (34) 

and Zi and I2 are the roots of the equation: 

+ (71 - 1 - to/2)Z + (2n2 - 4n — nm — m — 2) = (35) 

and Ai and A2 are two constants. 

Some special choices of n and m allows one to reconstruct the form of f{R) which are listed in Table®. 
Solution IV gives a negative sound velocity squared and thus is not acceptable. Solutions II, III, and V are 
acceptable for a thick shell. Again Israel junction conditions should be applied on it to have a complete 
solution. Finally solution I, with some values of n and m gives a physical acceptable solution for a star. 



3 Matching with exterior solution 

A look at Tables ([T|) and shows that only the solution I of Table (0) is appropriate for a star. To have 
a physical sphere of fluid surrounded by empty space, we set n < 2, so that the energy density and pressure 
decrease in a monotonic way with decreasing r, and the equation of state can be written as: 

p ~ Lop + Po{uj + I) (36) 

where uj — 2n+m-4: • Moreover in spite of obtaining infinite values for central density and pressure at r = 0, 
they are positive definite provided that po < 0, pi > and cj > 0. The pressure vanishes at the surface of the 
star r = Rs, and thus Rs — (-^)^^^"~^''. We have to have subluminal speed of sound and thus < w < 1 
and thus < m < 2 — n. 
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A{r) 



Table 2: The second class of solutions with X = X^r™ and F = For'" 
p{r) R{r) f{R) 



''s = tp 



4kXopi _ 
FqAo 

m(2w(w-l)-m(ri.+2)) 
n-2 



/o + xS?^'""^' forn^2 



m+2n— 4 



-4m^+4m-8n'^ + 16n 
-4m+4(n^-2n-l) 



/o + Fo-Rilni? for n = 2 



acceptable? 



some m, n 



II 



-314 I 30^1 
13r2 Xor* 



4(l-ofl+\/l+afl) /-l + Vl+aK 
3^^ V 1 

Where a = 30 (^)' 4^ 



5Xor^-52^i 
-4J(:or2 + 13Ai 



thick shell 



(See Discussion section) 



III 



-45 



30Ai 



SFqAi 1 



/o + 21n(-l + Vl + ai?)+ 
2v'l + ai? where a = 120 ( j^)^ 



Xor^ + 11.2Ai 
-Xor^-5.6Ai 



thick shell 



(See Discussion section) 



IV 



1GF„ 



-I 



V 



/InrS -I- 4,^4 „n - Mar2 I 12£iLli 1 

Aor +Air po iQ^r + inr 



^0 kXo 



-23 
5r2 



6Ai 



-ZXqt- +2aAi 
5Xor^-20Ai 



thick shell 



^r^-i^^lnr 



25^^gM ± I tanh-i ^1 + 6i?/a 



(See Discussion section) 



where a = (^) and 6 = 



24 A 1 
Xo 



Solution I specialize us to power law lagrangian plus a cosmological constant. Since f{R) models are of 
interest because of the fact that they can explain the late-time acceleration of the universe without the need 
for a cosmological constant, we set /o = 0. Fortunately, for this case in the absence of cosmological constant, 
there are two spherically symmetric vacuum static solutions in the literature. One is proposed in |12j and 
another is given by [5] . Therefore one can consider the matching of the interior metric to each of these forms 
of the exterior solutions at the surface of the sphere. Here we take the first choice. 

The exterior solution is thus[T^: 

Ae^ct = + Cr^ (37) 
Xe.t = 7^" (38) 

where C is a constant, a = ^"'"•^ , /3 = and 7 = 1— („ "2)2^" "^^ ■ The parameter n is subject 

to some constraints. The existence of chameleon-like mechanism (named by [33j so, to denote that the 
chameleon behavior in the Jordan frame differs slightly from the Einstein frame.), light deflection experiment 
and the perihelion precession of Mercury leads to the following bound on n [HI [T31 IMl IM] '■ 

- 1.6954 X 10"^^ <n< 6.7800 x 10"^'' (39) 

Using the exterior solution ((37|) and (p8| and the interior solution I of table ([2]), the junction conditions 
given by the equations (P5)) -(P7 )) . would be simplified into: 

A/iy =Q=^AoRT = Rs+ CR^, (40) 



Ai^: = =^ (m + 4) = 



(41) 



^0 

Ai? = =^ 7(5(m + 2)2 = (a + 2)2 (42) 

The first two equations determine the constants Aq and C in terms of i?s, n and m, while the third equation 
can be used to find m in terms of n. It can be solved exactly and there are two solutions: 

mi = \ > 43 

71 + 2 

and 

-2(n2-3n + 2) 

7712 = ^ -. 44 

71 — 4 

It has to be noted that only the ratio |^ is fixed according to the relation given in table ([2|). Therefore given 
a specific form of f{R) which fixes fo, one gets the constant pi. For 7711 one can simply see that pi = and 
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this wipes out the star solution. For m2 we find: 



^_4^o (n-l)3(n-4)4 

K (n-4)2(n2-2n-l) + (n-l)(n2-4)(n-3) ^ ' 



Also pq can be derived from the definition of Rs as: 

PO = -Pii?r' (46) 

As stated before, to have subluminal sound velocity we have to have < m < 2 — n. The acceptable 
values of n and m are plotted in Figure ([1]). Also taking into account the bound given by equation we 
see that values of n are very small and thus: 

5 

m ~ 1 n (47) 

4 

The last two junction conditions (j28p and (j29p determine the surface energy-momentum tensor. Equation 
(EH) gives: 

5 = (48) 

while equation (j28p simplifies to: 

5oo = A,F,K (^R7-' + 2^Rr-"y^ - ^V5^ (49) 

533 = 522 sin^ (51) 
5oo is the energy density and 522 and 533 are the tension densities on the star surface. 



4 Discussion 

In this paper we found some static spherically symmetric interior solutions of /(i?) theory of gravity. Of the 
eleven solutions, listed in table ([T]) and only one solution is physically acceptable for a star. In obtaining 
the solutions, we restricted ourselves to the solutions which have two properties. First, they are in the closed 
form and are obtained without any numerical or approximation method for solving the coupled high order 
differential equations of f{R) theory. Second, the corresponding form of f{R) function can be reconstructed. 
This means that we have not reported those solutions for which there are algebraic relations for R{r) and 
/(r), without the possibility for solving them to get f{R). 
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Figure 1: m as a function of n. Those values of m and n that hes in the shaded region lead to a subluminal 
sound velocity. 

Special attention is paid to the star solution, and it is matched to existing exterior solution and thus as a 
result one has a complete solution for a star and its exterior region for power law f{R) theories. The junction 
conditions arc applied on this star solution and the acceptable values of the constants are derived. 

It has to be noted that we are not able to apply the junction conditions for our thick shell solutions. This 
is because of the fact that we do not know the exterior solution for such /(i?)'s. Moreover those forms of 
f{R) are not of physical interest. 
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